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We introduce a mathematical framework for the analysis of the input-output dynamics of exter- 
nally driven memristors. We show that, under general assumptions, their dynamics comply with 
a Bernoulli differential equation and hence can be nonlinearly transformed into a formally solvable 
linear equation. The Bernoulli formalism, which applies to both charge- and flux-controlled mem- 
ristors when either current- or voltage-driven, can, in some cases, lead to expressions of the output 
of the device as an explicit function of the input. We apply our framework to obtain analytical 
solutions of the i — v characteristics of the recently proposed model of the Hewlett-Packard memris- 
tor under three different drives without the need for numerical simulations. Our explicit solutions 
allow us to identify a dimensionless lumped parameter that combines device-specific parameters 
with properties of the input drive. This parameter governs the memristive behavior of the device 
and, consequently, the amount of hysteresis in the i — v. We proceed further by defining formally a 
quantitative measure for the hysteresis of the device for which we obtain explicit formulas in terms 
of the aforementioned parameter and we discuss the applicability of the analysis for the design and 
analysis of memristor devices. 



I. INTRODUCTION 

According to classical electrical circuit theory, there 
are three fundamental passive circuit elements: the resis- 
tor, the inductor and the capacitor. Back in 1971, Chua 
challenged this established perception^. He realized that 
only five out of the six possible pairwise relations be- 
tween the four circuit variables (current, voltage, charge 
and magnetic flux) had been identified. Based on a sym- 
metry argument, he postulated mathematically the exis- 
tence of a fourth basic passive circuit element that would 
establish the missing link between charge and magnetic 
flux. The postulated element was named the memristor 
because, unlike a conventional ohmic element, its instan- 
taneous resistance depends on the entire history of the 
input (voltage or current) applied to the component 2 -! 3 -. 
Hence the memristor can be understood in simple terms 
as a non-linear resistor with memory. Chua and Kang^ 
generalized this concept to a broader family of non-linear 
dynamical systems, which they termed memristive, and 
demonstrated their usefulness in the modeling and un- 
derstanding of physical and biological systems such as 
discharge tubes, the thermistor, or the Hodgkin-Huxley 
circuit model of the neuron. Although in the interven- 
ing years experimental devices with characteristics sim- 
ilar to the memristor (i.e., hysteresis, zero crossing and 
memory) were investigated, researchers generally failed 
to associate or explain them in the context of memristive 
systems^. 

The memristor remained an elusive theoretical compo- 
nent until recently, when scientists at Hewlett-Packard 
(HP) fabricated a solid-state device which was recog- 
nized as a memristor— In that work, Williams, Strukov 
and co-workers also provided a simple heuristic model 
to understand how memristive behavior could emerge 



from the underlying physical process. The report of 
the fabrication of the HP memristor reignited the in- 
terest of the community in such circuits, due to their 
potential widespread applications. This has resulted in 
various attempts to build memristor devices based on 
different underlying physical principles, as well as ef- 
forts to understand pre-existing devices in the context of 
memristive systems. Examples include thin-films^^r— , 
nano-particle assemblies^, spintronics^ or neurobiolog- 
ical system a 16 i 17 . On the theoretical front, however, the 
analysis of these systems has been mostly restricted to 
numerics (i.e., temporal integration of model differential 
equations and numerical sweeping of parameters) due to 
the absence of a general mathematical framework that 
can provide analytical solutions to their dynamics. 

Here, we introduce a mathematical framework for the 
study of memristors that allows us to provide analytic 
solutions for their dynamics and input-output charac- 
teristics. The framework is based on the compliance 
of a general class of ideal memristor dynamics with Ja- 
cob Bernoulli's differential equation, a classic nonlinear 
equation that can be solved analytically. This formu- 
lation provides a powerful and systematic methodology 
for the analysis, characterization and design of devices 
governed by Bernoulli dynamics that does not rely on 
computationally expensive sweeping of parameters^—. 
Moreover, our analytical results can be used to reveal 
the relevant combinations of parameters that govern the 
response characteristics of the memristor. 

The paper is structured as follows. Firstly, we intro- 
duce the general theoretical framework for memristor dy- 
namics and the analytical solutions that follow from the 
Bernoulli differential equation. We then illustrate its ap- 
plication through the HP memristor—, which is shown 
to comply with the Bernoulli framework, and we obtain 
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analytical expressions of the output of the model when 
excited by three fundamental inputs: sinusoidal, bipolar 
square and triangular waves. To exploit further the in- 
sight provided by the framework, we define a quantitative 
measure of the hysteresis of the mcmristor in terms of the 
work done by the driving signal, and we use the derived 
analytical solutions to show that the hysteresis of the 
device depends on a specific dimcnsionlcss lumped quan- 
tity that combines all the parameters of the model. This 
lumped parameter relates fabrication and device proper- 
ties together with characteristics of the input signal and 
shows how such experimental parameters can be used to 
design and control the response of the memristor. 



II. MEMRISTOR DYNAMICS UNDER 
TIME- VARYING INPUTS 

A. Definitions and background 

Consider the four fundamental circuit variables: volt- 
age v, current i, charge q, and magnetic flux tp. There 
are six distinct relations linking these variables pairwise. 
Two of these relations correspond to the definitions of 
charge and magnetic flux as time-integrated variables: 



which serves us to define the memristance M. (<?, tpfa 



LI- 
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i(r)d,T 
v{r)dr. 



Three other links are given by the implicit equations that 
define the constitutive laws of the generalized fundamen- 
tal circuit elements: 

fn(v,i) = for the resistor, (f) 
fc(v,q) = for the capacitor, (2) 
fc(<f,i) = for the inductor. (3) 

In order to complete the symmetry of the system- 
theoretic structure, Chua's insight was to postulate that 
the remaining link between q and tp should be completed 
by another constitutive relation 



( 9f M \ 



dp 



i(t) = M(q,tp)i. 



(6) 



Since we are considering strictly passive elements, with 
M(q,<p) > 0, the constitutive law (j4]) defines a strictly 
monotonically increasing function in the (q, tp) plane. It 
then follows that we can express the memristor consti- 
tutive law as an explicit function either in terms of the 
flux, q = qjv[{v)i or i n terms of the charge, ip = tp M (q), 
indistinctly and to our convenience^!. 

Take, for instance, a current-driven mcmristor, in 
which i(t) is the time-dependent input under our control, 
then we can formally express the input-output relation- 
ship in terms of q, an integrated variable of the input i(t) 
over its past: 



v = M 



i{r)dT )i(t), 



(7) 



a form that makes explicit the memory of the device. 

Note, in contrast, that a time- invariant resistor gov- 
erned by ([U has a differential resistance TZ(v, i) 



dv 



w. 

\ dv 



di = lZ(v, i) di. 



(8) 



The input-output relationship of a current-driven passive 
resistor is then given by 



K(r})dr), 



(9) 



which is in general nonlinear but, unlike Eq. (J7J), it ex- 
hibits no memory. 

Hence, the memory property follows from the integra- 
tion necessary to go from the i — v plane, which is readily 
accessible to experiments, to the q — tp plane in which 
the memristor constitutive law is defined^. The mem- 
ristance (J5|) and the resistance © are thus dimensionally 
congruent but they only coincide if they are both con- 
stant, i.e., for the ohmic (linear and constant) resistor—. 



f M (v,q)=0, (4) 

which would correspond to a missing element: the mem- 
ristor. In this sense, the memristor complements the 
other three fundamental circuit elements as the fourth 
ideal passive two-terminal componentii££. 

Assume a time-invariant memristor with no explicit 
dependence on time, i.e., dfM/dt = 0. Then the time 
derivative of the constitutive relation ((4]) leads to 



dfM 


(df M " 


) ^ + 


(df M \ 


dt 


V dtp , 


K dt 


V dq J 



B. Solution of memristor dynamics under 
time-dependent inputs: the Bernoulli equation 

Although the constitutive relation of the memristor is 
defined in the q — tp plane, such devices are characterized 
experimentally in the i — v plane. Specifically, it is typi- 
cal to characterize their output dynamics in response to 
time-dependent inputs. The particular relationship im- 
posed by the definition of the memristor translates into 
a restricted form of the differential equations that gov- 
ern the dynamics in the i — v time domain. In particu- 
lar, the dynamics obey a Bernoulli differential equation 
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(BDE) that can be reduced to an associated linear time- 
dependent differential equation (LDE). In some cases, 
this structure can be exploited to provide analytical so- 
lutions for the input-output dynamics without the need 
for explicit numerical integration of the model. 
The Bernoulli equation has the forrr>22: 



and the resulting equation is: 



dx 



f(x)y = g(x)y r ' 



(10) 



where n £ N. If n — {0, 1}, the equation is linear. For 
n > 1 this family of nonlinear ODEs can be solved ana- 
lytically using the nonlinear change of variable z = 
to give the LDE 



dz 
dx 



+ (1 - n)f{x)z = (1 - n)g(x) 



from which the general solution follows: 
(1 — n) J m(x)g(x)dx + C 



II = 



m(x) 



(11) 



Here m{x) = e^ 1 ~ n - f( x ) dx is an integrating factor and 
C is the constant of integration. 

In the case of memristor dynamics, the mathematical 
transformation that relates the BDE/LDE pair is not 
only a convenient manipulation but is also associated 
with the experimental measurement setup and the in- 
trinsic description of the memristor: both the BDE and 
LDE can appear depending on the choice of externally- 
driven variable (i or v) in conjunction with the internal 
variable (q or p) through which the memristive behavior 
is controlled. 

In particular, consider an experiment in which the 
memristor is described in terms of the integrated input, 
e.g., a current-driven memristor with input i(t) whose 
memristance is charge-controlled, M(q). The dynamical 
system can then be written in the canonical form of a 
memristive system: 



q = i(t) 

v = M(q)i(t). 

This leads directly to the following LDE 
1 di(t)~\ \dM~ 



dv 
~di 



i(t) dt 



dq 



i(ty 



(12) 



(13) 



which has the standard analytical solution shown in Ta- 
ble u 

However, in certain situations, it might be more mean- 
ingful (or convenient) to express the memristance as a 
function of the integrated dependent variable. For in- 
stance, if the memristor is voltage-driven with input v(t) 
but the memristance is given in terms of the charge M (q), 
we have 



ip = v(t) 
M(q)i = v(t) 



di 


1 dv(i)~ 


i = — 


1 dM~ 


~dt ~ 


v(t) dt 




v(t) dq 



(15) 



which has a Bernoulli form. Indeed, this is the case for 
the first HP memristor model^, as shown in the next 
section. 

In general, the output dynamics of the general mem- 
ristor leads to four cases. Differentiating ^ with respect 
to time yields 



dv , .di 
— = M— + i 
dt dt 



dM . dM ' 

dq dp 



(16) 



If the device is current-driven, the dynamical equations 
are: 

dv _ I" 1 di(t)~\ \dM , dM 
~dt ~ 



i{t) dt 



dq 
dM 



M 



i{tf 



dp 

1 dM 1 



dq M 2 dp M 
while for the voltage-driven memristor we have: 



(17a) 
(17b) 



di 
dt 



1 dv{t) 
v(t) dt 



dM 1 dM I 



dq M 3 dp M? 



dM 1 dM 



dq M dp 



v{tf 
(18a) 

(18b) 



(14) 



where we have emphasized in the equations the ex- 
plicit dependence on time of the externally controlled 
drives. Note that the dynamics are formally either a BDE 
(Eqs. CE7p) and (fl8p)) or a LDE (Eqs. (fiTtt) and (pk)). 
which arc related through the nonlinear transformation 
given above. For strictly passive mcmristors^ (for which 
M(q), or equivalently M(p>), is always positive), the 
memristance can be represented as a function of either 
the independent or dependent variable, as convenient, 
due to the existence of the inverse function^. Therefore, 
once we specify the intrinsic integrated variable (i.e., q 
or p) that controls the specific memristor, we can obtain 
BDEs and LDEs as summarized in Table Q] with their 
corresponding general solutions. 

The above discussion shows that the constitutive rela- 
tionship of the memristor leads generally to Bernoulli dy- 
namics which can be reduced to the corresponding time- 
varying linear dynamics through a nonlinear transforma- 
tion. The analytical solutions obtained can be used to 
understand explicitly the parameter dependence of mem- 
ristor models as well as to aid in the design process with- 
out the need for parameter sweeping via numerical sim- 
ulations of the dynamics. In particular, we obtain below 
analytical expressions for the behavior of a memristor 
model when driven with different input waveforms. Fur- 
thermore, our analytical solutions allow us to show that 
the behavior of the system is fully characterized by a par- 
ticular lumped parameter that incorporates all the phys- 
ical parameters of the memristor. This renormalized pa- 
rameter is shown to encapsulate both the hysteretic and 
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Memristor Type 



Governing Differential Equation 



General Solution 



Current 
Driven 



Charge 
Controlled 



Flux 
Controlled 



Current 
Driven 



Voltage 
Driven 



LDE ^ 

dt 



1 di 
i dt 



dM{q) 



dq 



v(t) = i(t) 



Mr 



dM{g) 
dq 



(r)d7 



BDE ^- 

dt 


"1 di" 

i dt 


v — 


~dM{q) 1 
dq M 2 (q)_ 


v 2 v(t) = i(t) 


J dq M 2 (q) 


-i 



LDE * 
Voltage a* 
Driven j,- 

BDE - 
dt 



1 dw 
v dt 

1 

w dt 



dM{q) 1 



dg 7W 3 (g) 
dM(q) 1 



dg « 



w 2 i(t) = 



*(*) = «(*) 




dM(g) «(t) 
o dg 7W 3 (g) 

4 dX( g) 
dg 



dr 



-v(r)di 



LDE ^ 

dt 



1 di 
i dt 



dA4(^) 



dt/3 



M(<p) 



v{t) = i(t) 



dip 



BDE ^- 

dt 


"1 di' 

i dt 


v — 


'dA%) 1 
dip M{ip)_ 


v 2 v(t) = i(t) 




CdM{p) i(r) 1 
o dp M(tp) 


-i 


LDE 

dt 


T dv" 

v dt 


i = — 


'dM(p) 1 
dp M 2 (p)_ 


v 2 i(t) = v(t) 




f'dMip) v(r) d _ 
o dp M 2 {ip) 





di 

BDE — - 

dt 



1 efo 
w dt 



dX(^) 



d<p 



i(t) = 



^«(r)dr 



TABLE I. Governing differential equations and corresponding general solutions for different types of memristors, as characterized 
by their controlling variable, when driven by different inputs. 



nonlinear characteristics of the device and its memristivc 
character. 



C. Dynamics of the HP memristor under 
time-varying inputs 

We exemplify our approach with the HP model, which 
was introduced by Strukov et al. to describe the behav- 
ior of their fabricated memristor—. Their device consists 
of a thin-film semiconductor of TiC>2 with thickness D 
placed between two metal contacts made of Pt. The film 
has two regions: a doped region of thickness w with low 
resistivity Ro n due to the high concentration of dopants 
(positive oxygen ions) and an undoped region with thick- 
ness (D — w) and high resistivity Roff- The total re- 
sistance of the film is modeled as two variable resistors 
in series whose resistance depends on the position of the 
boundary between the doped and undoped regions. Ap- 
plying an external voltage across the device will have the 
effect of moving the boundary between the doped and 
undoped region effectively changing the total resistance 
of the device. Assuming ohmic electronic conductance 
and linear ionic drift with average vacancy mobility fi v , 
the device is modeled by the following pair of equations^: 



v(t) 
w(t) 



ON' 



D 



7Z-OFF 1 



D 



K-ON . u , 
Vv—p-l{t)- 



i(t) (19) 
(20) 



Equation (|20[) means that the position of the boundary 
between regions is proportional to the total charge that 
passes through the device. Hence Eq. (j!9[) is a memristor 



relationship of the form given in Eq. (|14[) with memris- 
tance controlled by charge^: 

kai \ -n ['R-on\ 2 TZoff\ 

M( q )=TZ FF + ^[-^) { 1 -K^)l 

= U FF-Kq- (21) 

When the device is voltage-driven (as is the case in the 
original model^) the dynamics is governed by Eq. (|15p . 
A simple inspection of Table fl] shows that the solution 
of the corresponding BDE leads to an i — v relationship 
with the following explicit form: 



v(t) 



^R 2 -2 K f*v(T)d7 



(22) 



where the constant of integration Rq € [Hon > Hoff] is 
the resistance of the device at t$, the initial time at which 
the input signal is applied. 

We can use the explicit solution (|22[) to calculate the 
output response of the HP memristor to different volt- 
age drives. Here, we will study periodic input signals 
v(t) with period To = 2-k/loq, amplitude A, and zero 

mean value, jr J Q v(t)cIt = 0. Such inputs are stan- 
dard test signals for electronic devices and lead to the 
classic characterization of the memristor through hys- 
teresis and Lissajous-type i — v curves^. In particular, 
we consider three test signals widely used in engineering 
settings, namely, the sinusoidal, bipolar piece- wise lin- 
ear and triangular waveforms, all defined over a period 
< t < T as follows: 



(i) Sinusoidal voltage input: 



a{t) = Asm(w t), 



(23) 
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(ii) Bipolar piece-wise linear wave: 



(ii) Bipolar piece-wise linear wave: 



n m (t) 



AJ_ 

m T 
A 

A / t 
m \T 

-A 

Ar_t 



o < 4 
— 1 1 



to < 



1) 



< to 

J£_ ^ 1 

To 

m < 



_/_ 

_/_ 

To 

m<4 



< 2 ~ m 



+ to < 



< g + TO 

< 1 — TO 

< 1 



(24) 



(iii) Triangular wave: 

a (f) =n m=1/4 (t). 



(25) 



For the bipolar piece- wise linear wave, the parameter < 
m < 1/4 determines the rise and fall time, assumed equal 
here. In particular, to = corresponds to the (Heaviside) 
square wave, while to = 1/4 is the triangular wave (|25[) . 
The three input waveforms arc shown in Figs. Ilajllcl 

It is convenient to normalize the i — v characteris- 
tics (l22l) in terms of rescaled variables: 



(26) 



a normalized time (or phase): x = ujq t 
a normalized input: v = v/A 

a normalized output: i = iRo/A, 



to give 



v(x) 



V 1 - I 3 Jo v( x ) dx 



(27) 



where we have defined /3, a positive dimensionless quan- 
tity which combines the parameters of the device and of 
the input drive: 



2A ( IZqn \ / T^-OFF 



r>2 I 



D 



\ T^ON 



-1 = 



2Ak 



(28) 



Note that the renormalized parameter /3 includes mate- 
rial and fabrication properties of the device (fi v , TZon, 
IZoff, D); parameters dependent on the preparation of 
the state of the device (i£o); as well as properties of the 
driving signal (A, ojq). We will show below that the mem- 
ristive behavior of the device is fully encapsulated by the 
parameter f3. 

Using the solution (j2"7) , the response to inputs (j2"3"|) - 
(|25|) over a period < x < 2ir is: 



(29) 



(i) Sinusoidal voltage input 
~ sin(x) 



— < 3- < 1 

y/l - /3(l - COs(x)) " 2n 



TCI 27T 



I = S 



V 1 


Bit { \ 




l 




- Png 2 {x) 


J_ 

m 


(1 _ JS.) 

\ 2 27T / 




- 




-1 






J_ 

771 





< < TO 

— Z7T 



2tt ^ 2 



^7T 



55 (x) 



i-TO<T^<i + TO 



i + to < tt- < 1 - m 



l_ m < x <1) 



(30) 



3i (x 
32(x 

33 (x 

34 (x 

35 (x 



(to — 2m 2 
(2-m-| 

f^-1) 2 



_x 1\ 

2tt 4 J 



(iii) Triangular wave: 
2 ii 



1 



2(1 -f) 



1-/3(2^-21-1) 
2(f -2) 

l- / 3(f-4a; + 47r) 



4 — 2-tt ^ 4 



(31) 



< 



< 1 



Figure [T] shows the response of the HP memristor for 
different values of f3 under the three periodic inputs 
shown in Figures [TalfTcl Figures lldlfTfl show the (non- 
linear) output current from the HP memristor for three 
values of /3, with the more positive values having a larger 
deviation from the linear output from an ohmic (linear) 
resistor with resistance Rq (dashed black lines). The re- 
sponse dynamics give rise to Lissajous-type i — v char- 
acteristics, which are typical of memristors when driven 
by a periodic excitation with zero mean valued. In such 
cases, the output is a double-valued function that forms 
a simple closed curve with no self intersections except at 
the origin. Equivalently, this corresponds to a hystere- 
sis loop crossing the origin in the i — v plane^ 5 -. As our 
results show, the amount of hysteresis in the device is 
dependent on /?, which changes as a combination of dif- 
ferent experimental factors. Figures |lg]llil illustrate the 
effect of j3 on the i — v characteristics: the hysteretic 
behavior of the device is enhanced for the most positive 
values of j3. Again, the dashed black line is the limiting 
case when the memristor tends to a linear resistor Rq, 
which corresponds to the limiting case of /3 = 0. 
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V(x) V(x) V(x) 

(g) (h) (i) 

FIG. 1. Input-output characteristics of the HP memristor under three inputs (arranged by columns) as given in Eqs. (1291) 
(|31[) , (a)-(c) The three input voltages as a function of the normalized time x: (a) sinusoidal input o(x); (b) bipolar piecewise 
linear wave n m=1 / g (a;); and (c) triangular wave A(x). (d)-(f) The effect of the parameter j3 on the output current of the HP 
memristor for the three inputs. In each figure, the normalized output of the device at different values of /3 (inset) is represented 
together with the normalized output of an ohmic (linear) resistor (black dashed line). As f3 approaches its upper bound /3 ma x, 
given in Eq. (|32[) , the output becomes more nonlinear and deviates from the linear response, (g)-(i) The effect of f3 on the 
i — v characteristics of the HP memristor for the three inputs. Each figure represents the plot of [v(x),i(x)] for different j3 
values together with the linear response of the ohmic resistor (black dashed line). The value of the linear resistor used here is 
Ro = 2f2. As /? — > /3 max , the i — v becomes more hysteretic. 



D. The lumped parameter /? encapsulates the 
dynamical response of the HP memristor 

Our analysis above shows that the dynamical response 
of the HP memristor is governed by the lumped parame- 
ter /?. This dimensionless parameter represents in a sin- 
gle quantity the collective effect of device parameters of 



diverse origin including: material properties, such as car- 
rier mobility and doping ratio (IZon, TZ-off, ^v)', device 
fabrication, such as the depth of the thin-film layer (D) ; 
preparation of the initial state of the device (Rq); as well 
as properties of the driving input (^4, ujq). Therefore, 
one does not need to consider each individual parameter 
separately — our analytical expressions show that similar 
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responses can be achieved by changing different proper- 
ties of the device if they lead to the same value of /?. 
Furthermore, for a given fabricated device one can tune 
the properties of the input drive in the experiment to 
produce particular output responses. 

We now examine some of the characteristics of the pa- 
rameter j3. ^From the definition (|2"5)) . it follows that 
P > since IZoff > TZon- As /? — > 0, the memristor 
tends to the behavior of a linear resistor with history- 
dependent resistance R . Hence, for a given fabricated 
device with particular physical characteristics, the lin- 
ear behavior will be revealed experimentally in the limit 
of low amplitude/high frequency drive (A — > and/or 
luq — > oo). In the process of designing memristors, lin- 
ear behavior will be more likely when the dimension of 
the device is large (D oo); when the mobility of the 
carriers is low (fi v —> 0); when the doping ratio is small 
(7Z-ON /TZ-off — > 1), or with a combination of all of those. 
In fact, note that j3 = when IZon = TZoff- This 
corresponds to building an ohmic resistor, rather than 
a memristor, with just one (undoped) region. In this 
case, the resistor always responds linearly with the same 
unique resistance independently of history or preparation 
protocol, i.e., Rq = IZoff- 

In the case of the HP model, the value of j3 is also 
bounded from above, as follows from requiring that the 
arguments in the square roots in Eqs. (j2"5)) - (f3"l"]) be pos- 
itive, so that the output current is real and finite. This 
leads to the following bounds for j3 for each type of input: 



< /3 < 1/2 for v(t) = cr(t) 

< $ < l/(?r - 2nm) for v(t) = n(t) 
< P < 2/ir for v(t) = A(i). 



(32) 



In order to make the comparison across input drives more 
direct, it is thus helpful to define the following rescalcd 
parameter: 



2/3 for v(t) = a(t) 

(it - 2Ttm)p for v(t) = l~l(f) 
(n/2)P for v(t) = A{t), 



(33) 



such that P € (0,1) for all three drives. When p = 0, 
the memristor becomes a linear resistor, while as P — > 
1, the HP memristor becomes maximally nonlinear and 
hysteretic and the separation of the two branches in the 
i — v plane is maximized. In the next section, we make 
this notion more precise through the introduction of a 
quantitative measure of hysteresis. 

It should be noted that the upper bound of P, that 
stems from the requirement that the output remains real 
and finite, is a manifestation of a limitation of the first 
simplified HP memristor model^, i.e., that the device 
does not saturate or break down as the memristance is 
driven towards its lower bound IZon- Although such a 
mechanism is not accounted for by Eqns. ([T9|) and (f20|). 
this limitation has already been addressed by models 
that contain a windowing function, in which it becomes 
increasingly difficult to change the memristance as the 



boundary between the doped and undoped region reaches 
either of the two ends of the device ^ 24 ' 25 . 



III. QUANTITATIVE MEASURE OF 
HYSTERESIS 



A. Definition 

Controlling and designing the hysteretic pinch of the 
i — v curve of the memristor is crucial for the use of these 
devices both individually or as part of larger circuits. It 
is important, therefore, to have a method to control the 
hysteretic effect of a particular model in order to achieve 
the required specifications for an application. We pro- 
pose here a quantitative measure of the hysteresis of the 
i — v curve in terms of the work carried out by the driv- 
ing input on the device. We will then use the analytical 
solutions obtained for the HP model in the previous sec- 
tion to calculate the hysteresis of the model for the three 
different drives in terms of the parameter p. Finally, we 
will show how these expressions may be used as an aid 
to fabricate a memristor with a prescribed i — v curve. 

Our measure of hysteresis is based on the calculation 
of the work done by the input signal through the device. 
Let H denote the positive scalar quantity measuring the 
difference between the work done while traversing the 
upper and lower branches of the hysteresis loop in the 
i — v plane: 

H = W + - W- = / ivdt- ivdt 
Jti Jt 2 

= / iddx- ivdx)= H, (34) 

Ro^o \J Xl Jx 2 ) Ro^a 

where t\ to ti and ti to £3 is the time required to move 
along the upper and lower branch of the hysteresis loop, 
respectively. Here, P(t) = i(t)v(t) is the instantaneous 
power and, as it is customary by now, H corresponds to 
the rescalcd quantity. 

Clearly, H is defined in terms of the energy dissipated 
by the device and becomes zero when the memristor 
tends to a linear resistor since W+ and W— coincide in 
that case. Therefore, it is useful to scale the hysteresis by 
Wo, the work done on the linear resistor Rq when driven 
for one complete cycle by the same signal for which H is 
evaluated: 



W * = TT I v 2 (r)dT = Ro f i 2 { 
Ro Jt Jt 



r)dr 



A 2 



v dx 



A 2 



Ro^o Ji-k Rqujq 
We define the scaled hysteresis as: 



Wo. 



H 



H 

Wo 



H 
Wo' 



(35) 
(36) 

(37) 
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B. Analytical expressions for the hysteresis of the 
HP memristor 



ar piecewise linear input 



We now use our analytical results for the HP model to 
obtain the hysteresis of this device under the three drives 
specified in the previous section. 



After some algebra, the scaled hysteresis of the bipolar 
input waveform v(t) = fl(t) is obtained as^ 7 -: 



1. Sinusoidal drive: 

Consider an HP memristor with sinusoidal input v(t) = 
a[t) and the corresponding output current i, given by 
Eq. (|29"]). The hysteresis H is then: 



H = 2 



sin 2 (x) 



yji - /3{1 - COS(X)) 



zdx 



(38) 



sin 2 (x) 



dx 



The work done by the input signal a(t) for a complete 
cycle on the linear resistor Rq is: 



H = 



3y/n? f 1 



0(3 - 8m) { p 



m arcsm 



i' In 



- 1 m' + f3(m — to') — a / /3to 



m'(l -/3) 




1-/3 




to' + (3(m — m') + y ml — j3m 



to'(1 - ^) - • 



\Jm' + P(jn-m') J L (41) 



Wo = / sin (x)dx = 7T. 
/o 



where to' = 1 — 2?n. 



Using the symmetry of the functions, the scaled hystere- 
sis H is given by: 



H = 



sin 2 (a;) 



zdx 



o y/1 - - cos(x)) 

sin 2 (a;) 



y/1 - j3(l - COS(X)) 



which can be rewritten a a 26 i 27 : 



£f 



32 y 7 !^ 



37T/3 2 



[2F(^,fc)-ii-(fc)] 



1-fj [E(k)-2E(<h,k)]-^\ 



1-/3 



(39) 



(40) 



Here, F(<j>,k) 



where fc 2 = S— and 02 = arcsin \ / — £ . 

P-l rz y 2-/3 

and £7(0, fc) are the incomplete elliptic integrals of the 
first and second kind, respectively, while A"(fc) and E(k) 
denote the complete elliptic integrals of the first and sec- 
ond kind, respectivel y 26 ^ 7 . 



P(x) 




FIG. 2. The instantaneous power, P — i v, as a function of 
the normalized time x for the HP memristor when driven by 
the bipolar square wave n m=1 / g (x). The output current is 
given by Eq. (|30[) . The areas w\ to wa indicate the work for 
the time intervals [0, 7r/2), [ir/2, n), [tv, 3ir/2) and [37r/2, 2ir], 
respectively. The hysteresis ()34[) is equivalent to H = W+ — 
W- = (u>2 + W3) — (wi + W4) where W+ and W- indicate the 
work along the upper and lower branches, respectively. The 
asymmetry in the instantaneous power is a reflection (and a 
measure) of the observed hysteresis. 
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H 




(a) (b) 



FIG. 3. (a) The normalized hysteresis H of the HP memristor model as a function of the memristor lumped parameter /3 under: 
a sinusoidal input (red solid line, Eq. (|40p ), a bipolar piecewise linear drive with m = 1/20 (green dotted line, Eq. ()4ip ) and a 
triangular input (blue dashed line, Eq. (|42p ~). The hysteresis curves achieve finite maxima /f ma x at j3 = 1. For a given lumped 
parameter (e.g., /3q — 0.9), we can obtain the hysteresis of the memristor when driven by any of the three inputs studied in 
Fig. [U The corresponding i — v characteristics of the HP model with /3o = 0.9 are shown in (b). Conversely, one can evaluate 
the hysteresis from i — v curves, obtained either theoretically or experimentally, and determine the corresponding /?. 



3. Triangular wave input 

The hysteresis for the triangular wave input v(t) = 
A(t) is obtained by particularizing Eq. (|4Tj) for m = 1/4: 



H 



1-/3 In 




(42) 



C. The dependence of the hysteresis on the 
parameter j3 

The definition of hysteresis (|34[) is based on the inte- 
gration of the instantaneous power consumed by the de- 
vice over the course of an input cycle. Figure [2] shows 
a plot of the instantaneous power for the HP model 
driven by the bipolar wave l~l(i). Each of the areas W\ 
to u>4 indicate the work done by the input signal for the 
corresponding time interval. The hysteresis is given by 
H = (w 2 + w 3 ) — (u>i + W4) where W + — w 2 + W3 and 
W- = w\ + W4. Note that u>2 = W3 ^ wi = W4. This 



asymmetry is a reflection of the difference in the work 
carried out on each of the two branches of the i — v curve 
due to the nonlinearity of the device and giving rise to 
the hysteresis. 

The expressions of the memristor hysteresis subject to 
the three drives given by Eqs. (|4U l) - (|42p are all explicit 
functions of /3, the lumped parameter that combines the 
physical parameters of the device together with the prop- 
erties of the drive. Figure |3a] shows that for all drives, the 
hysteresis of the device is zero when /3 = and increases 
as (3 — > 1. It is important to remark that the device has 
a finite maximum value of hysteresis it can exhibit, i.e., 
the value of H does not diverge as (3 — > 1. In fact, we 
can use our analytical expressions to calculate the upper 
bound of the hysteresis for each drive. For instance, from 
Eq. (|39[) the maximum hysteresis for the HP memristor 
driven by a sinusoidal input is: 



F max = H{p = l) = - 

IT 



1 



= -(1- 



sin 2 x 



cos(|) 
0.4972, 



dx — 



2 sin x 
cos(f) 



dx 



while for the triangular input, the upper bound of the 
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hysteresis is slightly lower: 

H0=l) = -j=(l- J) « 0.4552. 

Therefore, the maximum hysteresis exhibited by the HP 
memristor, understood as the difference between the 
work along the upper and lower branches of the i — v 
over a period of the input, is equivalent to ~ 50% of the 
energy dissipated by the equivalent ohmic resistor. 

The dependence of the hysteresis on the lumped pa- 
rameter P can be used during the design process of a 
memristor or during the characterization of a fabricated 
device. If the aim is to design a memristor with a pre- 
specified i — v response that needs to operate under a 
particular input, one may calculate first the hysteresis H 
of the desired i — v curve and the value of j3 that will 
produce the desired response. The identified P value can 
then be used to restrict our fabrication parameters. Sim- 
ilarly, for a given fabricated memristor, one can generate 
different i — v characteristics under a particular type of 
drive with varying frequency and/or amplitude, and/or 
under a different type of drive. For each i — v curve, 
the scaled hysteresis H can be obtained from the exper- 
imental data. If the data is well described by the HP 
model, our expressions could be used to fit some of the 
intrinsic parameters of the device that are contained in 
P. The interplay between hysteresis and j3 is exemplified 
in Fig. I3bl for the i — v curves generated by a memristor 
with Pq = 0.9 under the three input drives. 

IV. DISCUSSION 

In this paper, we have presented a mathematical frame- 
work for the analysis of the input-output dynamics of 
memristors. Although these are nonlinear elements, 
the form of the constitutive relation of the memristor 
leads in general to Bernoulli dynamics, which can be di- 
rectly linked to an associated linear differential equation 
through a nonlinear transformation. Table U highlights 
this general form of memristor dynamics and the duality 
that emerges when devices with distinct internal control 
variables arc interrogated with different externally con- 
trolled variables, i.e., charge- and flux-controlled mem- 
ristances which can be voltage- and current-driven. 

Our methodology allows us to obtain, in some cases, 
the output of the device as an explicit function of the 
input. We exemplified our analysis with the recently in- 
troduced HP memristor model^, which is shown to be 



a Bernoulli memristor, for which we obtain analytically 
the output current explicitly as a function of the voltage 
for three typical input signals, namely, the sinusoidal, 
bipolar piecewise linear and triangular waveforms. Our 
analysis led us to the identification of the dimensionless 
lumped parameter P, which combines physical parame- 
ters of the device as well as properties of the input drive, 
and controls the hysteretic properties of the i — v charac- 
teristics. We have shown elsewhere^ that the same pa- 
rameter also quantifies the amount of nonlinearity in the 
output spectrum of the device. Consequently, P controls 
both the hysteresis of the memristor and the harmonic 
distortion that the device introduces. The functional 
form (|27p makes apparent how p controls the mcmris- 
tive character of the device and reveals the fundamental 
interlinking between nonlinearity, hysteresis and memory 
in these devices. 

The explicit solutions thus obtained can be of use not 
only for the understanding of memristive dynamics but 
also as a means to parameterize experiments and to study 
the deviation from idealized models. They can also be 
used to devise experiments that can reveal the memris- 
tive properties of particular devices. For instance, it is 
possible to use P to design an input drive (e.g., its func- 
tional form, amplitude and frequency) that will enhance 
the hysteretic response of a particular memristor, depen- 
dent on the physical properties of the device and the un- 
derlying transport mechanism of the charge carriers. The 
parameter P also indicates which fabrication parameters 
are likely to enhance (or reduce) memristive behavior. 
In the HP model, the small dimensions of the device, a 
large carrier mobility, and large differences between the 
doped and undoped resistances all contribute to make the 
memristive behavior more conspicuous. The results pre- 
sented could be applicable to experiments given that the 
simple HP model studied here has already been shown 
to provide a reasonable approximation of the behavior of 
particular experimental realizations^. 

Although wc have focused here on particular exam- 
ples, this analytical approach has been extended to other 
memristor models^. Furthermore, equivalent analytical 
input-output relations can be obtained for networks of 
memristors and other mem-elements connected to circuit 
parasitica^. The study of such directions as well as the 
investigation of more realistic memristor models (in par- 
ticular those that account for non-linear dopant kinetics) 
will be the object of further research. 
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